Electronic Interface Reconstruction at Polar-Nonpolar Mott Insulator Heterojunctions 
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We report on a theoretical study of the electronic interface reconstruction (EIR) induced by po- 
larity discontinuity at a heterojunction between a polar and a nonpolar Mott insulators, and of the 
two-dimensional strongly-correlated electron systems (2DSCESs) which accompany the reconstruc- 
tion. We derive an expression for the minimum number of polar layers required to drive the EIR, 
and discuss key parameters of the heterojunction system which control 2DSCES properties. The 
role of strong correlations in enhancing confinement at the interface is emphasized. 

PACS numbers: 72.80.Ga,73.20.-r,71.10.Fd 



I. INTRODUCTION 

Two-dimensional (2D) electron systems have been a 
fertile source of interesting physics over the past few 
decades, playing host to the fractional and integer quan- 
tum Hall effects and cuprate superconductivity among 
other phenomena. The most widely studied and most 
thoroughly understood 2D electron systems are those 
that occur near semiconductor heterojunctions. In these 
systems carrier densities and disorder strengths can be 
adjusted using modulation doping and the electric field 
effect, and high sample quality can be achieved using lat- 
tice matched materials and epitaxial growth techniques. 
These 2D systems are well described by Fermi liquid the- 
ory, at least at magnetic field _B = 0. Rapid recent 
progress in the epitaxial growth of complex transition- 
metal oxidesi^i^i^i^i^ii»2, foreshadows the birth of an en- 
tirely new class of 2D electron systems, one in which 
electronic correlations are strong even at i? = and non- 
Fermi-liquid behavior is common. 

The present authors have recently argued^ that 
the modulation doping technique commonly used in 
semiconductor heterojunction systems can also be ap- 
plied to complex-oxide heterojunctions to create high- 
quality two-dimensional strongly-correlated electron sys- 
tems (2DSCESs) localized near the interface between 
two different Mott insulators. Our study was based 
on a generalized Hubbard model which captures some 
key features of these systems, and on combined insights 
achieved by applying Hartree-Fock theory (HFT), dy- 
namical mean-field theory (DMFT) and Thomas-Fermi 
theory (TFT) approaches. This theoretical strategy 
echoed that adopted by Okamoto and Millisi^ in address- 
ing the two-dimensional electron systems which can occur 
near the interface between a Mott insulator and a band 
insulator. In both cases the properties of the 2D electron 
systems reflect a delicate balance between space-charge 
electric fields and strong local correlations. 

These 2D systems are clearly unusual in many respects. 
For instance, their quasiparticles interact with the 3D 
spin-excitations of the Mott insulator barrier materials. 
Since researchers now appear to be at the cusp of achiev- 
ing experimental control over this kind of 2D electron 
system, it is interesting to explore the possibilities theo- 



retically in an attempt to find useful theoretical frame- 
works, and hopefully also to partially anticipate some of 
the new phenomena likely to occur. To mention one ex- 
ample, it has recently been suggested that a spin-liquid 
insulator could be found in a Mott insulators bilayer— 
system with an appropriate doping profile. 

The apparent electronic interface reconstruction (EIR) 
discovered^i^ at the heterojunction between the band 
insulators LaAlOa and SrTiOa is unique in complex ox- 
ide interface studies and does not have an analog in semi- 
conductor heterojunction systems, although there is still 
some uncertainty about the role of oxygen vacancies in 
the measured conductivity. For example, it has been 
proposed^^ that the oxygen vacancies are responsible for 
the insulating behavior observed at the p-type interface. 
The reconstruction is forced by the polarity difference 
between LaAlOa and SrTiOa. Polarity discontinuities 
are normally weaker at semiconductor heterojunctions 
and, when present, the electric fields to which they give 
rise are usually screened by relaxation of atoms near the 
interfacei^. In complex transition metal oxides, however, 
it was discovered^ili, that the polarity discontinuity field 
can be screened electronically by transferring electrons 
between surface and interface layers, changing the va- 
lence of transition metal ions in these layers. Electrons 
dodge the polar catastrophe^ without essential atomic as- 
sistance. It has been showrt^ that such an interface can 
be tuned between metallic and insulating states by rela- 
tively modest electric fields, suggesting the possibility of 
novel electronic devices. 

In this paper we examine EIR at a heterostructure be- 
tween model polar and nonpolar Mott insulators. We 
have in mind for these effects Mott insulators with cu- 
bic peroskite structures. Possible materials combinations 
are discussed in Section VI. The long-ranged Coulom- 
bic space-charge fields are treated within the Hartree ap- 
proximation, and strong on-site Coulomb interactions are 
treated by solving HFT, DMFT, and TFT equations. We 
find that a minimum number of polar layers is required 
to induce an EIR because of strong local correlations, 
whereas the number of nonpolar layers does not play 
an essential role in determining the electron distribution. 
We also find that when a 2D electron system is present 
near the heterojunction, it is strongly localized because 
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FIG. 1: Schematic illustration of our model for a thin film 
containing a polar-nonpolar Mott insulator heterostructure. 
The model consists of Np layers of a polar Mott insulator 
perovskite {AMO^) and A*' — A'^p layers of a nonpolar Mott 
insulator perovskite {A' M' O3.). The z direction is chosen to 
be the layer-by-layer growth direction and the charge density 
is assumed to be uniform in each x — y plane. The symbol A 
represents a group III element with nominal valence -1-3 in the 
AO layer while the symbol A' represents a group II element 
with nominal valence +2 in the A'O plane. The AO layer 
therefore has a surface charge density of -1-1 per A atom while 
the A'O layer is neutral. The symbols M and M' represent 
transition metal ions in group IV and V respectively. Both M 
and A'/' have id} configuration before the electronic interface 
reconstruction (EIR). 



of strong correlations which limit the potential of the 
electronic system to lower its energy by spreading out in 
the growth direction. TFT results provide a simple way 
of understanding the more microscopic DMFT and HFT 
calculations and are remarkably accurate. 



II. SINGLE BAND HUBBARD MODEL 

Fig. [1] illustrates the model heterostructure we inves- 
tigate in this paper. We consider a thin film composed 
of Np polar perovskite Mott insulator layers {AMOz) 
and N — Np nonpolar perovskite Mott insulator layers 
{A'M'03), where A {A') is a group III (II) element with 
a valence of +3 (+2) in the AO{A'0') layer. M (M') is 
a group IV {V) element which is a Mott insulator 
because of strong local repulsion among the d-orbitals. 
In this study we disregard the interesting complications 



associated with orbital degeneracy and do not account di- 
rectly for hybridization between the transition metal and 
oxygen ions. We therefore use a single band Hubbard 
model to describe the d valence electrons. Long-ranged 
Coulomb interactions are described realistically by ac- 
counting for charges at the oxygen sites and on the A, 
A', and M, and M' ions. The resulting model Hamilto- 
nian is: 
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where i labels a metal site and (grouping coplanar oxy- 
gen and A cation charges) = 1 in the present case. 
The d-orbital hopping term is expressed in Eq.([T]) as a 
sum over links; most of our calculations have been per- 
formed for a model with nearest-neighbor hopping. Sec- 
ond nearest-neighbor hopping has little influence on our 
main results, but will impact the 2DSCES Fermi surfaces 
as discussed below. The sizes and alignments of the gaps 
in the two Mott insulators are fixed by the interaction 
and site energy parameters: U{zi) = Up and ed{zi) = Cdi 
for the polar layers and U{zi) — Unp and ed{zi) = 6^2 
for the nonpolar layers. In an ideal cubic perovskite 
unit cell with lattice constant a, Ri = a{ni,mi, Zi) and 
Rf'"^' = a{n, + 1/2, + 1/2, Zi + 1/2) respectively. The 
third term in Hcoui accounts for a net effective charge 
per M' atom in the M'02 plane; this charge is the sum 
of the oxygen ion charges and the +5 charge of the M' 
ion when in a d^ configuration. In this work we have 
neglected the uneven distribution of charge within each 
layer. In this calculation we for the most part assume 
that Cd is constant, but emphasize that the difference 
between d-orbital energies on opposite sides of the het- 
erojunction can be an important parameter. As in case 
of semiconductors, band lineups play a key role in de- 
termining hetero junction properties and are likely to be 
difficult to estimate with sufficient precision. Accurate 
values will most likely have to be determined experimen- 
tally for each heterojunction system considered. The re- 
quirement of overall electrical neutrality in the thin film 
implies that the d orbital occupancy, averaged over the 
thin film p, satisfies p = 1 . 

Throughout this study we use the Hartree approxi- 
mation for the long-ranged Coulomb interaction and ne- 
glect the corresponding exchange term, since our goals in 
studying this toy model are purely qualitative. The self- 
consistently determined Hartree mean-field Hamiltonian 
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where pj — '^^{n-ja) is the electron density on site j, 
/ = for jz < Np, and / = 1 for jz > Np. In the 
continuum Umit, H'^J^i becomes: 
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III. HARTREE-FOCK THEORY 

In HFT the Hubbard U term is replaced by its mean- 
field counterpart: rii^nii ^^^("'i-o") ?iio-- As pointed 
out^iii^ previously, there are almost always a number of 
distinct self-consistent solutions of the HF equations, cor- 
responding to various local-minima and saddle-points of 
the HP energy functional. Here we present results for 
solutions with G-type antiferromagnetic (AFM) order, 
or (tt, tt, tt) spin order, and ferromagnetic (FM) order. 
These solutions are representative of ordered states which 
occur frequently in complex transition metal oxides and 
are often global minima of the HF energy functional. The 
philosophy of following this procedure is that HFT can 
neither reliably judge the competition between different 
types of order, nor predict whether or not long-range or- 
der succumbs to fluctuations. The results nevertheless 
very often provides a reasonable description of local cor- 
relations in the strongly correlated regime. 

Our HFT results for Np = 2 and Np = 5 are summa- 
rized in Fig. [5] For the calculations reported here, the 
total layer number = 10 and typical values were cho- 
sen for the interaction parameters. {Up/t — UNp/t = 20 
and Uc/t = 0.8.)ia For both AFM and FM states FIR 
occurs only for the model thin film with the larger value 
of Np. This feature is a result of the competition between 
the polar catastrophe^, and strong local correlations. Al- 
though the polarity discontinuity at the interface favors 
a transfer of electrons from the outmost surface to the 
interface, this charge rearrangement is opposed by the 
Mott-Hubbard gap. It is instructive to consider space- 
charge induced band bending diagrams like the cartoon 
illustration of Fig. [3] which contrasts the Np = 2 and 
Np = 5 cases. For Np = 2, the bending of the upper and 
lower Hubbard bands is not large enough to force the 
chemical potential p at the heterojunction outside the 
gap region. When more polar layers are grown, fi even- 
tually crosses the lower Hubbard band of the outermost 
polar layer and the upper Hubbard band of the nonpolar 
layer nearest the interface, forcing the occurrence of an 
EIR. If the top layer of the polar material was a posi- 
tively charged AO layer rather than a negatively charged 
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FIG. 2: Electron charge distributions for AFM and FM 
Hartree-Fock states with Up/t = Unp/t = 20, Uc/t = 0.8, 
and ed constant for (a) Np — 2 and (b) A'p = 5 cases. The to- 
tal number of layers is TV = 10. EIR occurs only in the model 
thin film with the larger number of polar layers. Confinement 
of the 2DSCES is stronger in the AFM state solution because 
of its stronger local correlations. 



MO2 layer, the charge transfer would occur from the 
lower Hubbard band near the heterojunction to the up- 
per Hubbard band near the surface. 

An analogous competition has already been dis- 
cussed and studied experimentally in LaAlOa/SrTiOs 
heterostructures^. In that case the energy which com- 
petes with the polar catastrophe is the energy difference 
between the top of oxygen p-bands in nonpolar SrTiOs 
and the bottom of aluminum d bands in LaAlOa, in- 
stead of the correlation gap in the example studied here. 
Since it always costs the gap energy to add more elec- 
trons on any layer of the heterostructure, electrons accu- 
mulate initially near the interface to gain the most elec- 
trostatic energy, eventually spreading out through the 
heterostructure to minimize kinetic energy cost. Quan- 
tum confinement of the 2DSCES associated with EIR is 
stronger for states with stronger on-site correlations be- 
cause the kinetic energy is then a smaller component of 
the overall energy budget. This is the reason that in Fig. 
[21[b) electrons are more confined to the interface in the 
AFM state than in the FM state. (The FM state has 
larger band widths and smaller correlation gaps, at least 
within HFT.) Based on the discussion above, we expect 
that for fixed Up and Unp, the minimum number of po- 
lar layers required to achieve reconstruction iV™ will be 
smaller when Uc is larger. N^ will also tend to decrease 
when the Hubbard gap is reduced in either layer; in Fig. 
|3]we illustrate the dependence on Unp when Up is fixed. 
Note that a smaller value of Unp favors ferromagnetic 
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FIG. 3: (Color online) Cartoon illustration of band bend- 
ing for (a) Np = 2 (b) Np = 5. The support of the upper 
Hubbard band's spectral weight is indicated by yellow shad- 
ing while the support of the lower Hubbard band's spectral 
weight is indicated by solid red. These illustrations reflect 
the spectral weight only at the M and M' transition metal 
sites of M{M')02 layers with the end point on the polar side 
of the interface at the last MO2 layer and the starting point 
on the nonpolar side of the interface at the first M'02 layer. 
Since the last MO2 and the first M'02 layers are separated 
by a AO layer, the electrostatic potential reaches its mini- 
mum at the first M'02 layer. (The electrostatic potential is 
of course continuous along the growth direction, as required 
by the Poisson equation.) (a) ^ lies in the gap region and 
no EIR occurs, (b) /i crosses the lower Hubbard band on the 
polar side and upper Hubbard band on the nonpolar side, in- 
dicating the appearance of EIR. If the top layer of the polar 
material was a positively charged AO layer rather than a neg- 
atively charged MO2 layer, the charge transfer would occur 
from the lower Hubbard band near the hetero junction to the 
upper Hubbard band near the surface. 



HF states. 

A more microscopic view of EIR is provided in Fig [5] 
where we illustrate the electronic structure of the AFM 
solution of the Hartree-Fock equations for a thin film 
with Np = 5. The parameters used to construct this 
illustration are the same as those used for the Np = 5 
case in FigH^b). The iV = 10 thin film has 20 two-fold 
spin-degenerate 2D bands because of the reduced trans- 
lational symmetry of the AFM state, and charge neu- 
trality requires that 10 bands be occupied. In Fig. El^a) 
we can identify 10 relatively widely spaced bands which 
correspond at high energies to the upper Hubbard band 
and at low energies to the lower Hubbard band in the 
polar material. The bands are widely spaced because 
of the average electric field in the polar material; for 
the same reason these bands are quite highly localized 
in individual atomic layers. The two groups of more nar- 
rowly spaced bands correspond respectively to the upper 
Hubbard band and the lower Hubbard band of the non- 
polar material. Because the space-charge electric fields 
are much either almost fully or at least partially screened 
out in the non-polar material by charge accumulation at 
the interface, the bands are closely spaced and the corre- 
sponding eigenstates contain more inter-layer characters. 
When EIR occurs there is a weakly avoided anticross- 
ing between the top-most lower Hubbard band state in 
the polar material and the upper Hubbard band states 
of the non-polar material. This anticrossing is apparent 
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FIG. 4: Ground-state phase diagram (a) vs. Np and Uc/t for 
Up = Unp = 20 and (b) vs. Np and Unp for Up = 20 and 
Uc/t = 0.8. The total number of layers is fixed at = 10. 
(a) For fixed Up and Unp, EIR occurs for a smaller number 
of polar layers if Uc is larger. Note that the FM states always 
have higher ground-state energy than AFM states in these 
parameter regions, (b) For fixed Up and Uc, EIR requires 
more number of polar layers when Unp is larger. As indicated 
here, a smaller Unp also favors FM states over AFM states in 
Hartree-Fock theory. For both (a) and (b) the phase boundary 
for the occurrence of EIR is accurately reproduced by the 
Thomas-Fermi approximation, Eq[7] 

in FiglSI^b) which expands the region of the electronic 
structure close to the Fermi energy. In the ground state 
the polar lower-Hubbard band state is lower in energy 
at the Brillouin-zone center but higher in energy toward 
the Brillouin-zone edges. (For a model with only nearest- 
neighbor hopping the lower upper band has its maximum 
and the upper Hubbard band has its minimum along the 
line from (tt, 0) to (0,7r).) 

For the parameters used in this calculation (with sec- 
ond nearest-neighbor hopping t' /t ~ —0.15) the EIR 
leads to two partially occupied 2D bands illustrated in 
Fig. [DJc), a nearly full lower-Hubbard band very local- 
ized on the surface and with two inequivalent hole pock- 
ets in the AFM Brillouin zone centered on (7r/2, it/2) and 
(7r/2, — 7r/2) and a nearly empty band concentrated on 
the first non-polar layer with an electron pocket centered 
on (77,0). 



IV. THOMAS-FERMI THEORY 

Layered oxide materials have a natural Thomas-Fermi 
approximation in which the total energy is expressed as 
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FIG. 5: (Color online) Electronic structure of the A'p = 
5 AFM state whose charge distribution was illustrated in 
FiglJ^b). This plot was calculated with a non-zero second 
nearest-neighbor hopping parameter t' /t = —0.15. Energies 
are measured from the chemical potential of the thin film, 
(a) All 20 2D bands along high symmetry lines of the AFM 
Brillouin zone, (b) Crossing between the polar material sur- 
face layer lower Hubbard band and the five non-polar mate- 
rial upper Hubbard band layers. The space charge physics 
implies that these six bands will be close to the chemical po- 
tential. Note that only two bands are partially occupied, (c) 
Dependence of the total band occupancy on position in the 
Brillouin-zone on a color scale with occupancy 18 in red, 20 
in green, and 22 and in blue. 



the sum of the total electrostatic energy and a local- 
density-approximation for the band and correlation en- 
ergies. For layered structures the contribution to the en- 
ergy from each atomic layer is approximated by the en- 
ergy per layer of an electrically neutral 3D system with 
the average site occupancy of that layer. Minimizing this 
energy with a fixed total density constraint leads to the 
following TP equation: 



Y^\z-zj,\-Y^{p{z')-f)\z'-z\ (5) 



and /ig is the overall chemical potential. Following the 
same approach described in Ref^, we can solve the TF 
equations for different ordered states using calcu- 
lated from HF solutions for the three-dimensional Hub- 
bard model. We confirmed that TFT very accurately re- 
produces the results of HFT for the charge distribution 
between the layers, as in the case of modulation doping-. 

One advantage of using this TF equation is that we 
can derive some key properties analytically. For example, 
if we assume that the d-orbitals of polar and non-polar 
sides are both described by one-band Hubbard models 
and that the charge transfer occurs from a single surface 
layer to a single interface layer Eq.(|l]) reduces to 

Edi +A'p(l ^ '5"-) + 27rC/cA^p(l - 2 5n) = €^2+ ^J-Np{^ + Sn) 

(6) 

where the Hubbard model chemical potentials are ex- 
pressed as function of density in units of electrons per 
atom, Sn is the charge transfer, and the on-site energies 
which are usually set to zero in Hubbard models, restored 
on each side of the heterojunction. Since the left hand 
side is a monotonically decreasing function of Sn and the 
right hand side is a monotonically increasing function of 
dn, this equation has a solution and EIR will occur if 

2TrUcNp > [ed2 + /ijvp(l+)] - [edi + A'p(l-)] = Eig (7) 

i.e. if the electric potential drop across the polar layer ex- 
ceeds the spatially indirect gap Ejg between the bottom 
of the upper Hubbard band in the non-polar material and 
the top of the lower Hubbard band in the polar material. 
Since the variation of chemical potential with density is 
strongly reduced outside of the gap in strongly correlated 
material, it follows that the amount of charge transfered 
once UcNp exceeds the minimum value is given accurately 
by 



2t:UcNp 
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Eig 



(8) 



(4) 



As the right hand side of Eq. [8] becomes comparable to 
0.5 at larger values of A^p, the fundamental assumption of 
our simple model (namely that the polar catastrophe is 
avoided by electronic rather than atomic reconstruction) 
becomes increasingly suspect. When EIR occurs it is 
likely to lead to relatively small 2D Fermi surface pockets, 
as illustrated in Fig. [5] 



V. DYNAMICAL MEAN-FIELD THEORY 

Hartree-Fock theory can describe the physics of strong 
local correlations only when the correlations are static 
and lead to EM or AFM ordered states; HF solu- 
tions for PM states completely miss local correlation 
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FIG. 6: Electron distribution calculated by HFT and DMFT 
for PM state with (a) Np = 2 and (b) Np = 5. Other param- 
eters are the same as those used in Fig [2] The comparison 
between HFT and DMFT for PM state demonstrates the role 
of on-site correlation in the EIR. The solutions of DMFT show 
that EIR only occurs with more polar layers and the confine- 
ment is stronger compared to HF solutions. 

and do not reproduce the Mott-Hubbard gap. DMFT 
captures correctly the insulating gap created by on- 
site correlations even for inhomogeneous multilayered 
systems^iiSiii^. For this reason DMFT calculations for 
paramagnetic states (PM) provide important additional 
insight into EIR in Mott-Hubbard insulator heterojunc- 
tion systems. For the impurity solver necessary in appli- 
cations of DMFT, we adopt the simple two-site method 
proposed by M. Potthoff^^, which has been successfully 
applied to a variety of different transition metal oxide 
heterostructures^iiS,. The results for Np = 2 and Np = 5 
are presented in Fig. [6] with Up/t = UNp/t = 20 and 
Uc/t — 0.8. When a paramagnetic state is described 
by HFT, charge transfer from the surface to the inter- 
face occurs (incorrectly) even for Np = 2. For Np = 5, 
the screening of the surface layer electric field incorrectly 
begins immediately in the polar material. The DMFT 
calculations demonstrate that a larger value of Np is re- 
quired even if the spin-degree of freedom is not ordered 
in the Mott insulators. These results confirm that a min- 
imum thickness of polar layers is required for EIR when 
on-site correlations are strong. 



VI. MATERIAL CONSIDERATIONS 

It is not so easy to determine on purely theoretical 
grounds which material combinations could be used to 
realize the EIR proposed here. In the family of cubic 



perovskites, candidate d polar materials include LaTiOa 
and YTiOa and possible non-polar materials include 
SrVOa include CaVOa. The band offsets of the vari- 
ous possible materials combinations Cdi — td2 can be es- 
timated from existing electronic structure calculations. 
LDA calculations based on Wannicr functionsi^ indicate 
that the t2g bands tend to be higher in energy relative 
to the oxygen p levels in LaTiOa and YTiOa than in 
CaVOs and SrVOa. Assumming that the oxygen p lev- 
els are closely aligned across the heterostructure, we can 
expect that LaTiOa and YTiOa will have higher d or- 
bital energies than CaVOs and SrVOs. Another interest- 
ing observation from the LDA calculations is that energy 
differences between different t2g orbital associated with 
distortions of the ideal cubic perovskite structure tend to 
increase across the series SrVOa-CaVOs-LaTiOs-YTiOs. 
These additional features related to the t2g orbital de- 
gree of freedom may bring in new physics which can not 
be explored in our single-band thin film model. More- 
over, the tendency toward magnetic ordering may also 
be influenced by the presence or absence of t2g orbital 
degeneracy, which will in turn be influenced by reduced 
symmetry near the heterojunction. The same calcula- 
tions show that CaVOa and SrVOa are metallic in the 
bulk. Lower coordination and electric fields that reduce 
inter-layer hopping are likely to tip the balance toward 
the insulating state in thin films and especially near het- 
erojunctions, so these materials are still candidates for 
realization of the physics studied here. 

Although sophisticated first principle calculations can 
provide some insight and will help with the construction 
of realistic phenomenological models, the consequences of 
orbital degeneracy or near-degeneracy for complex oxide 
heterojunctions may depend on subtle issues of many- 
particle physics. Another aspect not captured realisti- 
cally in this qualitative study is the dielectric constant 
which in general could be a very complicated function of 
positions throughout the heterostructurcs^^. Experimen- 
tal information when available will likely play an essential 
role in achieving a full understanding. 

One of the most important effect of the band lineups 
in our model is its influence on where the electrons ac- 
cumulate when EIR occurs. We have shown that in the 
case of constant electron accumulation near the het- 
erojunction appears on the first M'02 layer. This will 
change if for the polar materials is low enough so that 
the sum of ea and electrostatic potential is minimized at 
the last MO2 layer, instead of the first Af O2 layer. 



VII. SUMMARY 

In this paper we have studied electronic interface re- 
construction (EIR) in polar-nonpolar Mott insulator het- 
erostructure systems. We use a simple one-band Hub- 
bard model supplemented with long-range Coulomb in- 
teractions which assumes that the polar catastrophe in 
layered polar materials is avoided by electronic rather 
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than lattice reconstruction. We find that EIR occurs 
only if a thin film contains a minimum number of po- 
lar layers N^; when model paramaters representative of 
perovskite oxides arc chosen N™ ~ 4. In this paper, 
the electronic structure of polar-non-polar Mott insula- 
tor heterojunctions was investigated using Hartree-Fock 
theory (HFT), Thomas- Fermi theory (TFT), and dynam- 
ical mean- field-theory (DMFT). Using TFT we are able 
to derive an analytical expression for the dependence of 
EIR on model parameters, which may be useful in se- 
lecting suitable materials. For polar and non-polar Hub- 
bard interaction parameters {Up and Unp respectively) 
a stronger Coulomb interaction parameter Uc will reduce 
iV™, while for Up and Uc fixed, larger Unp will increase 
iV™. Strong on-site correlations enhance quantum con- 
finement of the interface states which accompany EIR, 
resulting in 2D systems that are nearly completely con- 
fined to an atomic plane. When the Mott insulator bar- 
rier materials have antiferromagnetic order, Hartree-Fock 
calculations demonstrate the emergence to 2D electron 
systems with Fermi surface areas and characters that de- 
pend sensitively on the details of the model. The quasi- 
particles of this system will differ from ordinary 2D elec- 
tron systems because they will interact strongly with the 
3D spin excitations of the antiferromagnetic Mott insu- 
lators. The unusual properties suggest the likelihood of 
interesting non- Fermi liquid states in some regions of the 
model's parameter space. 

We have also discussed some possible material combi- 



nations which might lead to the EIR we envision and the 
key role of the alignment of d-orbital energies across the 
heterojunction. Some realistic features of perovskites, 
for example, orbital degeneracy and lattice distortions, 
are not captured by the simplified model used in this 
work. We note that the exchange term from the long- 
ranged Coulomb interaction, which we have neglected in 
this study, should enhance the correlation gap. As a re- 
sult, the charge distribution may be somewhat different if 
this term is also treated self-consistently. These features 
of realistic systems are likely to play an important role in 
practice, especially so when the charge transfer from the 
surface to the interface region is large enough to sensibly 
alter the chemical bonding as discussed previously for the 
Mott-Insulator/Band-Insulator heterojunction casein. 

In our view, it will likely never be practical to attempt 
a complete and accurate description of complex oxide 
interfaces of any type from first principles. Insights from 
ab inito theory, simplified tight-binding models, and most 
critically from experiment arc all likely to play a role in 
achieving robust, predictive, phenomenological models of 
these systems. This work is intented as a step toward 
establish the correct form of useful phenomenologies. 
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